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Abstract — We present a unified method to compute the 
exact performance of ARQ schemes for wireless networks. 
A wireless channel is modeled here as a stochastic sequential 
machine. Accordingly, we propose a channel state ordered 
representation of the performance provided by an underly- 
ing block code. Based on proper events introduced to assess 
reliability, throughput, and delay performance of an ARQ 
scheme, we employ a discrete system model to obtain a ma- 
trix description for ARQ schemes used on wireless channels. 
Then, closed formulae are developed for all examined per- 
formance characteristics. Numerical results are given for a 
wireless ATM system. 

I. Introduction 

Wireless networks can be characterized as either wire- 
less local area networks mainly designed for operating in an 
indoor environment [1][2J or wide area networks [3]. The 
asynchronous transfer mode (ATM) is expected to play an 
important role in future generation wireless networks as well 
[4]. Wireless networks are supposed to support high qual- 
ity voice, video, and high speed data services. Therefore, 
one of the challenges when developing a wireless network for 
highly demanding applications is to combat the complex er- 
ror mechanisms of the radio path to guarantee high quality 
services. To support the services provided by a wireless net- 
work, sophisticated error control techniques have to be used. 

In order to maintain adequate system performance, block 
codes are used in combination with a retransmission strat- 
egy establishing an automatic repeat request (ARQ) scheme. 
The fundamental ARQ schemes are classified into stop-and- 
wait (SW), go-back- N (GBN), and selective repeat (SR) 
ARQ. Clearly, a performance analysis of an ARQ scheme 
has to include a suitable channel model. 

Fading channels such as in wireless networks cause nonin- 
dependent channel errors, so that a memoryless channel is 
not a realistic model. Because of statistical dependencies be- 
tween errors, a performance analysis of an ARQ scheme in a 
wireless channel requires more computational effort as in the 
case of a simple memoryless channel. Contributions dealing 
with ARQ schemes in nonindependent errors are very limited 
and focus on specific types of channel errors [5]. Assump- 
tions taken to simplify the mathematical tractability often 
result in a model close to a memoryless setting, e.g. sym- 
metrically dependent errors [6]. 

This paper presents an analytical method for evaluat- 
ing performance of ARQ schemes in nonindependent errors, 
i.e. channels with memory. An algebraically decoded lin- 
ear block code is utilized for error detection. Correspond- 



ing to the used stochastic finite-state channel model, per- 
formance of the employed block code is structured into a 
channel state ordered representation as well. By defining 
suitable events to be assessed, using automata theory for dis- 
crete system modeling and the concept of matrix functions, 
the proposed matrix approach leads to closed formulae for 
reliability, throughput, and delay. The unified method is ap- 
plicable for any generative discrete channel model and gives 
exact results. 

The paper is organized as follows. In Section II, we present 
a discrete system model for ARQ schemes in wireless chan- 
nels. On this basis, closed formulae for reliability, through- 
put, and delay performance are developed in Section III, IV, 
and V, respectively. Numerical results are given in Section 
VI and Section VII concludes the paper. 

II. Modeling of ARQ Schemes in 
Wireless Channels 

In this paper, a wireless channel is modeled as a stochastic 
sequential machine (SSM). We make use of the stationary 
probability distribution 



on the S channel states 5 and the state transition matrix 

D = [p(*'W]sxs, (2) 

where p(s'\s) denotes the probability when the channel is in 
state s that the subsequent state will be s'. 

A block code is utilized for error detection and is alge- 
braically decoded. Usually, performance of a block coded 
system without any retransmission strategy is characterized 
in terms of 



Pc 
Pe 
Pd 



with 



Probability of correct decoding, 
Probability of erroneous decoding, 
Probability of error detection, 



P c + P c +P d = 1. 



(3) 



Performance of a block code can be further classified by the 
probability P D of decoding and the probability P f of an 
error. We have 



Po = Pc + P e and P f = P e + Pd> 



(4) 
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To assess the performance of an ARQ scheme in a wireless 
channel, we assume the performance of a block code is given 
in a channel state ordered representation, e.g. Pc(s^\s^) } 
where s*°* denotes an initial channel state and s (n) is the 
channel state reached after transmitting a packet of length 
n. This requirement can be satisfied by using the method 
presented in [7]. The various conditional probabilities are 
taken as elements of matrix probabilities 



Pc = [Pc(« ( "V 0) )]sxS, 
P, = [p«(s (n) |s (0) )]5xS, 
Pd = tp d (s (,, V 0) )]5xS, 

which are interconnected via [8] 

P c + P e + P d = D". 



(5) 



(6) 



A matrix equivalent to (4) can be formally written as 

P D = P C + P e = [pd( S ( "V 0) )]sxs, (7) 



P f = P e + P d 



[Pf(s {n) \sW)) SxS . 



(8) 



With this approach, it is possible to include the chan- 
nel state transitions during the round-trip delay by means 
of the state transition matrix. In doing so, the elementary 
events to be considered are the decoding outcomes after the 
v ih transmission of a packet and the number of state transi- 
tions occurring between consecutive retransmissions of that 
packet. We use the notation: 

c v : correct decoding after the v ih transmission, 

e v : erroneous decoding after the v th transmission, 

d v : error detection after the v th transmission, (9) 

sJJ 1 : m channel state transitions between the v ih and 
{v + l) th transmission, 

where m represents round-trip delay in bits. Let 5 be the 
data rate, r the round-trip delay time, n the code or packet 
length, and let the operator \x] denote the smallest integer 
equal or greater than x. Then, the parameter m is given by 



for SW, 

for SR and GBN. 



(10) 



Hereinafter, we assume a symmetric forward channel, an 
error-free feedback channel, and perfect synchronization. As 
a consequence, an analysis of the impact of the channel error 
process on the all-zero codeword is sufficient and neither 
acknowledgment nor frame losses occur. 

III. Reliability Analysis 

Let £ be the event that a decoder releases an erroneous 
packet, either after the first or after some retransmissions of 
that packet. We introduce a subevent 

£ v : Decoder releases an erroneous packet after it 
has been transmitted v times. 



Each subevent £ v is an intersection {a} n {&} = {a, b] of 
elementary events (9) and given by 

£v={d u s?, d 2i 4 n , cUl C-i, (H) 
Thus, £ may be rewritten as union of subevents £ u as 



£ = {j£u 



(12) 



Therefore, reliability of an ARQ scheme can be classified 
indirectly by way of the probability P{£) of an undetected 
error, i.e. the probability of occurrence of the event £. The 
higher the probability P(£) is, the lower is the reliability 
and vice versa. 

Using the concept of matrix probabilities, the stochastic 
nature of a subevent £ v can be specified analytically by 



P(£ v ) = (PrfDT^.Pe. 



(13) 



Because the subevents in (12) are pairwise disjoint, the 
matrix probability of an undetected error is simply given by 
the sum 



P(f ) = £P(£„) = ^CPdD" 1 )- 1 • P e 



(14) 



Since the matrix probability Po of decoding and the ma- 
trix probability P c of correct decoding can be easily com- 
puted by using the method proposed in [7], it is convenient 
to substitute P<* and P e in (14) by 



P d = D" - P D and P e = P 0 - P c . 



(15) 



Performing some operations on (14), the matrix probabil- 
ity of undetected error of an ARQ scheme in a channel with 
memory can be expressed in closed form as 



P(£) = (I-D n+m +P 2) D TO )- 1 (P D -P c ) 



(16) 



Finally, we sum over all channel states and as a result 
obtain the desired scalar probability of undetected error 



P(£) = * 0 -P{£)e, 



(17) 



where e denotes an all-one column vector. 

Note that for an infinite round- trip delay and under the 
assumption of a regular Markov chain of channel states, each 
row of the state transition matrix is given by the stationary 
probability distribution [8]: 



(18) 



Accordingly, for an infinite round-trip delay equation (17) 
can be rewritten as 



Pe,SSM 



Pc,SSM + Pe.SSM * 



(19) 



where P c ,ssm and Pe,ssM denote the probability of correct 
and erroneous decoding of a block code used on a SSM, 
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respectively. In other words, statistical dependencies be- 
tween errors within a packet are still included in the anal- 
ysis whereas correlation between subsequent transmissions 
vanishes. 

Similarly, for the special case of a binary symmetric chan- 
nel (BSC) equation (17) approaches the known result [9]: 



Pe,BSC 
Pc,BSC + Pe.BSC ' 



(20) 



where P c ,dsc and P e ,BSC denote the probability of correct 
and erroneous decoding of a code on a BSC, respectively. 

IV. Throughput Analysis 

Throughput of any ARQ scheme is dependent on the num- 
ber of packets which could be transmitted during the average 
time spent to decode a particular packet for the first time. 
Hence, we focus on the random variable 

T: Number of retransmissions, including the first 
transmission, until a packet is decoded either 
correctly or incorrectly for the first time. 



A particular outcome v of T will be denoted as 



(21) 



and can be written as an intersection of elementary events, 
that is 

T v = R,sr><*2,sJ\ — , , » D v = c„ U e„}. (22) 

In order to compute the average number of packets which 
could be transmitted until a packet is decoded for the first 
time, we introduce an expected matrix 



mr = E{T} = [E{T}(*< n >|*<°>)] 



SxS 



(23) 



where the traditional concept of an expected value is gener- 
alized and structured into channel states. On that basis, we 
first examine throughput in channel state ordered represen- 
tation and then average the result over all states: 



mr = (To * mr * e — > 77 ~ m T 



(24) 



For that purpose, we need the matrix probabilities related 
to the outcomes T Ui v = 1, • * ■ , 00. These are given by 



P{%) = (P d D m ) v - 1 -P D . 



(25) 



Because for SW ARQ and SR ARQ only packets detected 
in error have to be retransmitted but for GBN ARQ also a 
finite number of subsequent packets must be retransmitted, 
we have to distinguish between 



Tarq(%) 



for SW and SR, 
for GBN, 



(26) 



with subscript ARQ e {SW, GBN, SR}, and where the 
parameter N= \(S • r)/n] -1 represents round-trip delay in 
packets. Then, we have 

00 

^„ = £ T ARQ {%) ■ (PrfDT-'Pa. (27) 

By manipulating (27), the expected values for the number 
of transmissions until a packet is decoded for the first time 
can be obtained in closed form as 

m rsw = to * F^ 2 P D • e, (28) 

m TGBN = <To-[F» l +NF- D 2 (D n -P D )D m ]P D .e, (29) 

mTsr = ^0 ■ F^ 2 Pd • e, (30) 

where the following matrix has been introduced to simplify 
the notation: 



F D = I-D n+m + P D D" 



(31) 



Finally, the desired throughput of an SW ARQ, GBN 
ARQ, and SR ARQ scheme used on a channel with memory 
are given in closed form by 



VSW = 


(n + S • r) • m Tsw 




k 


r}GBN = 










k 


VSR = 






» " m T S R ' 



(32) 
(33) 
(34) 



V. Delay Analysis 

The delay associated with an ARQ scheme can be de- 
composed into queueing delay, transmission delay, and rese- 
quencing delay. In this paper, we focus on the transmission 
delay as result of the impact of a stochastic error process on 
data transmission. First, we characterize delay by the ex- 
pected value of the number of transmissions until a packet 
is decoded correctly for the first time and the associated 
standard deviation. Then, we can simply include round-trip 
delay time and the processing time spent at the receiving 
end for decoding. 

For that purpose, we introduce a random variable 

V : Number of retransmissions, including the first 
transmission, until a packet is decoded cor- 
rectly for the first time. 



A particular outcome v of V will be denoted by 
V v = {V = 1/}. 



(35) 



We examine delay and delay variation by counting only 
correctly decoded packets. We assume that a genie assists 
the decoding process to request a retransmission for a de- 
coded packet in the case it is released in error by the de- 
coder. 
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In this context, a particular outcome of V can be written 
as intersection of elementar y events (9) and is given by 

V» = {h, s?, / 2 , s 2 m > Z^C-i, <v}, (36) 

where {/„} = {e„} U {d v }. Thus, the matrix probability of 
a particular outcome becomes 

P(2>„) = (P/DT 1 -^ (37) 

Hence, the average number of retransmissions until a 
packet is decoded correctly for the first time is character- 
ized by the expected matrix 

CO 

mp = E{X>} = ^i/-(P / D m )"- 1 .P c . (38) 

For the sake of a set of homogeneous performance formu- 
lae, we substitute 

P, = D n -P c (39) 

in (38) and by further examining this sum, one may straight- 
forwardly show 

m v = cr 0 (I-D n+m +P c D m )- 2 P c .e. (40) 

Similarly, we are going to represent delay variation by the 
standard deviation 

a v = >jE{V*)-ml. (41) 

Because the average number of retransmissions can .be 
evaluated by (40), we only need to compute the moment 
matrix 

CO 

E(D 2 ) = • (P/D' 71 )"- 1 • P c . (42) 

As a result, we finally obtain the required moment 

E{V*) = a 0 -(2F~ 3 -F c - 2 )P c .e, (43) 

where 

F c = I-D n+m + P c D™ (44) 

Having developed closed formulae for the average num- 
ber of retransmissions m-u and the standard deviation av, 
we can now include the physical dimension time into the 
analysis. We are going to classify delay time r into three 
fractions: 

T = Tbiock + Trtd + Tdec, (45) 

where 

mock = n/<$ : Delay time caused by organizing sequences 

of bits into packets, 
Trtd ~m/6 : Round-trip delay time, 
Tdec ' Processing time spent for decoding. 

On that basis, we simply have to multiply the average num- 
ber of retransmissions and the associated standard deviation 
by the delay time and obtain the average delay time and de- 
lay variation: 

tu t ~ r ■ m-o and o T = r • ux>- (46) 



VI. Numerical Results 

We examine a GBN ARQ scheme used on a two-state 
Gilbert-Elliott channel (GEC). The GEC is defined by the 
transition probability P = Prob(B\G) and Q = Prob(G\B) 
between a 'good' state G and a 'bad' state B, respectively, as 
well as the bit error rates pa and p% in the respective state. 
An effective BSC is taken for comparison having the same 
average bit error rate p b as the GEC. In the following, the 
concept of a slow fading channel will be used to emphasize 
when transitions between states are rather unlikely. The 
corresponding GEC parameters are chosen as 

P = HT 6 , Q = 3 - lO" 6 , p G = lO" 6 , p b = f(p B ). 

Suppose the number of information bits taken is k = 424, 
which may reflect the number of bits to be protected in a 
wireless ATM system. A shortened cyclic code is employed 
for error detection and defined by the generator polynomial 

g(x) = x 16 +x 15 +x 13 +x 9 +x 7 + x 6 +x 5 +x 3 +x+l. 

Note that the round-trip delay m is used as parameter in 
all presented plots. A round-trip delay in a range up to some 
100 bits may apply for an indoor wireless network whereas 
higher round- trip delay corresponds to mobile radio systems. 

Figure 1 shows the probability P(£) of undetected error 
of the (440,424) block coded GBN ARQ scheme versus av- 
erage bit error rate. The scheme is used on a slow fading 
channel. Obviously, P(£) approaches to the probability P e 
of erroneous decoding of the underlying block code (dashed 
line) when we increase the round- trip delay. In other words, 
neglect of statistical dependencies between channel errors in 
adjacent packets would result in an overestimated reliability. 

Figure 2 compares the throughput of the examined GBN 
ARQ on a slow fading channel and a BSC. When transi- 
tions between channel states are very unlikely, the decrease 
in throughput in high average bit error rates turns out to be 
more severe on the GEC compared to the performance on 
the effective BSC. As far as throughput is concerned, the su- 
perior reliability in high average bit error rates on the GEC 
leads to no improvement. This is because in the interval of 
high average bit error rates throughput has already dropped 
to almost zero. 

We eventually consider average delay time versus average 
bit error rate. We assume a data rate of 5 = 2Mbps, which 
gives T bfock = 220/is and r rtd = 0 to 0.55 for m = 0 bit to 
10° bit. The delay fraction r dec has been neglected. Figure 3 
shows the results for the average delay time. As it can be 
seen from the plot, for average bit error rates p b < 10~ 3 , 
the GBN ARQ provides the same constant delay on both 
channels. On the contrary, at high bit error rates the de- 
lay increases exponentially on the effective BSC but changes 
only to a higher niveau of constant average delay time for the 
GEC. The reason is that the probability of correct decoding 
on an effective BSC approaches zero for high average bit er- 
ror rates. On the GEC the probability of correct decoding 
can be maintained non-zero because of the contribution in 
the good state G. Especially for the slow fading channel, 
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the stationary state probabilities for being in the states G 
and B are 75% and 25%, respectively. In addition, the av- 
erage periods of time spent in G and B are 10 6 and 1/3 • 10 6 
discrete time cycles, respectively. 

VII. Conclusion 

We have presented an analytical method to evaluate the 
exact performance of ARQ schemes in wireless channels. 
Corresponding to the employed stochastic finite-state chan- 
nel model, performance of a block code is structured into 
a channel state ordered representation as well. Events as- 
sociated with reliability, throughput, and delay of an ARQ 
scheme in a wireless channel have been introduced and are 
taken as a foundation for an exact performance evaluation. 
By using a discrete system model in combination with ma- 
trix functions, closed formulae are obtained for all examined 
performance characteristics. Because of the employed ma- 
trix approach, an implementation of the proposed method 
on a digital computer is straightforward. 
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Figure 1. P(S) of a (440,424) GBN ARQ scheme 
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Figure 2. 7] of a (440,424) GBN ARQ scheme 
Parameter: Round-trip Delay ra 
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